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Abstract 

It is shown that high accelerating gradient can be obtained in a specially constructed 

system of electron (positron) bunches, moving in cold plasma with definite density. These 

combined bunch systems do not generate the wake fields behind them and can pass 

through the plasma column in a periodic sequence. The consideration is carried out 

numericaly and analyticaly in one dimensional approach, (which can be applied to finite 

system when its transverse dimensions are larger than plasma wave length, divided by 
• »— i . 

27r). The possibilities of the experimental tests by measuring the predicted energy gain 



> 



are discussed on the examples of Argonne Wakefield Accelerator and induction linac with 
typical parameters. 

1 INTRODUCTION 

Langmuir already noticed, that the beam, which has passed the plasma column, contains 
a significant portion of electrons with the energies higher than the initial energy of the 
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beam. In the recent times various groups (see e.g. 0-0) also observed experimentally 
the effect of the selfacceleration. 

The possibility to accelerate the electrons of the bunches, moving in plasma, was 
considered numericaly in |3J, || and some attempts of the analytical treatment of the 
problem have been performed in ||- ||. 

In the work [lQj|, which is based on the previous results |L1| used as a zero order ap- 
proximation, nonlinear dynamics of the one- dimensional non rigid ultrarelativistic bunch 
of electrons, moving in the cold plasma, is considered by multiple scales perturbative 
approach. Square root of the inverse Lorentz factor of the bunch electrons is taken as a 
small parameter. 

It was shown [ID], in particular, that for specialy constructed system consisted from 
two bunches, first one dense enough ^> ^n Q , and the second one with low density 
n£ <C |no, the electrons on the rare side of the second bunch, where electric field E 
is negative and large, can be significantly accelerated. In the first approximation of 
multiple scales method, the change of the momentum of the bunches electrons is given 
by Apb = —eEt, where t is a acceleration (deacceleration) time interval which, due to 
applicability of the multiple scales method, is t < uj~ 1 ^ 1 / 2 . 

The essential restriction on the acceleration rate in the considered cases comes from 



the steady state assumption adopted in [IDJ. If plasma electrons momenta in the wake 
field behind the first dense relativistic bunch, moving in underdense plasma, exciedes same 
critical value, allowed by steady state assumption, the plasma electron density behind 
the bunch turn out to be negative, which probably means that wake wave breaking limit 
is achieved. Restriction on the plasma electrons momenta inside the bunch automatically 
restricts the maximum value of the generated electric field. 

In the present work another proper combination of the charged particles (e 1 * 1 ) bunches 
is proposed, which is constructed in such a way, that it does not generate the wake 
fields behind the combined bunch at all. The plasma electrons momenta, inside the 
combined bunch being negative, remains in the allowed domain and their maximum 
value is restricted only by the total momentum conservation low. The maximum of the 
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subsequent electric field then also can be large enough. 

The combined bunch in the considered case consists (at least) of three parts:first one 
is a relativistic negative particles bunch, with density nj^ ^> |no, which generates in a 
plasma inside the bunch strong electric field, which brakes the bunch particles (generator). 
The second bunch is a low density bunch with particles njf 1 < |no (it can be negative 
particles bunch, positive particles bunch), which inverts the sign of the electric field 
generated by the first bunch (invertor). The third bunch with the parameters, which can 
coincide with that of the first one, provided on its rear part the zero values of electric 
field and plasma electrons momenta, coincided with the same values on the front of 
the first bunch (damper). It means that plasma behind the system of bunches remains 
unperturbed and no wake field is generated. 

Acceleration of the particles of the rear part of the second and of the third bunches, 
where electric field is negative, will take place. The goal of the following calculations is 
to estimate the maximum possible values of the acceleration parameter. 

In the next section of the work the vivid approach to the problem, based on the nu- 
merical solutions of the subsequent nonlinear equations and their graphic representations 
on phase plane are presented. The third section devoted to the analytical treatment of 
the problem, which based on the exact solutions of the equations and provides the general 
expressions for the energy gain, acceleration rate and acceleration length, the lengths of 
the different parts of the combined bunch as a functions of charge densities and Lorentz- 
factor of the bunches. Finally, the numerical examples are presented, which consider the 
possibilities of the experimental observations of the predicted selfacceleration of charged 
particles, in particular, on Argonne Wakefield Accelerator (AWA) and induction linacs. 

2 PHASE PLANE IMAGES OF THE COMBINED BUNCHES, MOVING IN COLD PLASMA 

Consider a plane rigid electron bunch with the infinite transverse dimensions and length 
[0, d] along the z-axis. Velocity of the bunch in the lab system is v z = v o ((3 — ^), 
plasma is cold with the immobile ions. The charge density of the plasma electrons is 
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n^the case, when n& > |no called underdense regime, rib < |no-overdense regime. The 
plasma electrons motion equations described in hydrodynamic approximation is 

W + ^ = ~E,or ( A - A) )_ = -E; (1) 

Continuity eq. for plasma electrons is 

|^ + = 0, or A[ n ' e (/3 e -/5)] = 0, (2) 

and Maxwell (Coulomb) equation for the electric field is 

dE' E' 



dz dz' 

In eqs. (§0) 



i i i 
n n -n- n h 



(3) 



,/ . / , /4vre 2 n b 
t = cjpt z = k p z, Up = I 



rn 



mavp , , n e . rib 



E = (47mmc 2 ) 1/2 £' = -^^E',n' e = -,n' t 



> >% — — 

e n n 



Pe = p —^ e = v — , 7 = (i-/5 2 )- 1/2 ; 

mc c 

5 = 2 — f ot (steady state approximation, n- is an arbitrary density, which finaly will be 
chosen) . 

Continuity eq. (0) has an integral 

n e {(3-(3 e )=n (3- (4) 

constant is defined from the boundary condition, z = d, (3 e = 0, n e = n®. In ([|) and in 
what follows the superscripts "prime" are omitted. 
From eq. (Hj) it follows that 

£ = = = ^T? e -P Pe (5) 

and eqs (@)-(@) can be written in the form 

*''-^ V-%^ =a ~ 72 (6) 

(see also 0, H -[£!). 



Eq. (^|) has an "energy" integral 

S = i$ ,2 + 7 2 [$ - /3($ 2 - 7~ 2 )] -a$ (7) 

where a = n^/r^the boundary conditions at z = d are p e = 0, $ = 1,$' = and 
gives £ = 1 — a, but in what follows the eq. (|7|) considered at the arbitrary boundary 
conditions. 

The integral (0) allowed to interpret equations as equations for the point with 
unit mass with the "coordinate" $ and "velocity" $' moving in the potential U = Uq + 
U x = 7 2 [$ - ($ 2 - 7~ 2 ) 1/2 ] - a$. 

On Fig. 1 the {7 as a function of $ is drawn for different values of a and for 
7 = 10. The negative values of a = correspond to the bunch of positive charged 
particles, positive a corresponds to the bunch of negative charged particles. For a < 
0, 7~ x < $ < 1 and $ > lfor a > 0, when boundary conditions, when z = d, are 
$ = 1 and $' = 0. From Fig. 1 it is seen that for a < the motion is always 
periodic, it is also periodic for < a < 1/2, and non periodic for a > 1/2 JO]] . 
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Fig. 3. Phase-plane portrait for a=0.4, 7=10. Separatrix: 2=9.96 . Fig. 4. Phase-plane portrait for a=10, j=10. Separatrix: 2=9. 

Fig. 2, 3, 4 represent the phase plane (<£>', <£>) portrait of the system of cold plasma-rigid 
charged particle bunch in the cases, when a = —10, 0.4, 10 and 7 = 10 at different 
boundary conditions (different £). The closed phase trajectories correspond to periodic 
motion. 

Trajectories on Fig. 2, 3, 4 can serve as a building blocks to construct the phase 
portraits of the combined bunches, moving in cold plasma. It must be represented by a 
closed loop, starts from the boundary condition $ = 1, $' = and accomplished at the 
same point $ = 1, $' = on the phase plane. An example of such a loop is presented on 
Fig. 5. The loop on Fig. 5b from the point : $ = 1, = (p e = 0, v e = 0, E = 0, 
when z = d) and the curve OA represents the motion of the bunch of negative charged 
particles (generator) = 10, when at the rear side of the bunch very strong but positive 
electric field E = — <&' can be generated. Taking this point as a new boundary conditions 
for the second bunch of positive charged particles (invertor), the trajectory is turned to 
the point B, where electric field is zero and then to the point A' where the field —E is 
obtained. From point A' to the point the trajectory described the motion of the third 
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bunch of the negatively charged particles (damper) on the rear side of which the initial 
boundary conditions exist. It means that no wake field excites and the plasma behind the 
third bunch remains unperturbed. As a consequence, the above mentioned restriction on 
the value of the field E , which arouse from the condition of the stability of the wake field, 
is removed. 
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Fig. 5a. Uniform ee + e' bunch. Patameteres oi the bunch ate: 
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Fig. 5b. Phase-plane portraitfor e'e'e' bunch. 



This is not only the condition for the applicability of steady state regime for the 
description of the zero order approximation in the multiple scales approach, but also 
open the possibility to repeat more than once the procedure, sending the next combined 
bunches into the unperturbed by previous bunches plasma. 

The Lorentz factor of all three bunches in the presented case is, of course, the same 
7 = 100, and 



n 



(2) 



Jl_ = _i y 3 > = = ^_ = io. 

n n 



(3) 



electron bunches are ~ 15.77^ each. 

Z7T 



The total length of the positron bunch is 28.57^ and lengths of the first and third 

On Fig. 5a the potential $ and electric field E — — 1| as a function of z are presented. 
It is evident from Fig. 5a that positrons of the first part of the second bunch, where E > 0, 
can be accelerated as well as all electrons of the third bunch (E < 0). 

The similar picture may be obtained by construction of the combination of three 
bunches with negatively charged particles (electrons). First bunch is dense nj^ /n 3> 1, 



second has a low density \n Q — n b 3> 7 2 , and third is again dense n h 3> n Q . 




dj=2999.65*y2n. Fig. 6b. Phase-plane portrait for combined electron bunch. 

The example of such a combination is presented on Fig. 6a, 6b, where electric field and 
phase portrait are drawn for the case, when 

n W = n f ] = 10n ,n ( b 2) = 0.4n ,7 = 100. 

It is seen from Fig. 6a, that the length d,2 of the second bunch-invertor is too large 
in considered case. Acceleration of the bunches electrons take place when E < i.e. 
mostly in the rear side of the second bunch-invertor and, at some extend, in the third 
bunch-damper. 

The role of the bunch-invertor may play also the plasma itself (n b =0). In such a 
case the electric field interacts on plasma electrons, which are moving in the direction 
opposite to that of bunch, with the initial velocity equal to zero. So they will accelerate 
in the region where E > 0, and drag in the region where E < 0. The electrons of the 
third bunch, where E > will accelerate. The charge distribution inside the bunches can 
be nonuniform. 

On Figs. 7a, 7b the results of the numerical solution of eqs. are displayed for 

the case of sinusoidal charge distribution inside the electron- positron-electron bunches 
span by plasma in-between. The bunches have the equal lengths d = 15 7p, devided by 
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plasma columns with lengths do = 60^,7 = 100, 




It is neccessary to note that the wake field behind the combined bunch due to small in- 
accuracies, presented in actual calculations practically exists, but the amplitude of this 
wake is small and far from wave breacking limit. This last condition defines the tolerances 
on the parameters of combined bunch and must be taken into account in actual calcula- 
tions for the proposed experimental tests. In the case of the sequence of the combined 
bunches, passing through plasma, the conditions, which define the tolerances on bunch 
parameters, will be more severe. 

3 ANALYTICAL APPROACH 

In this section the results of the exact analytical solutions of eqs. for the com- 

bined bunch with the uniform charge distributions inside the constituting bunches are 
presented. Results for single rigid (electron) bunch with the zero boundary conditions 
p e = 0, E = 0, when z — d) are obtained in ||11|| . 

Boundary conditions at the end of the first electron bunch moving in underdense 
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plasma rib k n o are p e = p , E = E , where 



E = ± [2 (n?> - n )] 1/2 {(1 + p 2 ) 1 / 2 - 1 - a^} 1 ' 2 



a 



(i) 



'/. — '"oyj ) • 1 ■ - /'o ' '/'o f (8) 

nl /n\ - n 

The field inside the second electron bunch, moving in overdense regime njf* < l/2n 
with the same velocity is 

£ (2) = ± [2 (n - n?>)] V2 {a - (1 + p 2 ) - , 



a = 



,(2) 



no— n 

(l) (2) 



T2T 



< 1, a< 2 )/3 < 1, 



0) 



When p - *■ a — > 1, when 



p < 0, | po| > 1 



n« - nf 2 > 



n - n 6 (2) 

When plasma electron momenta is equal to 



(l + /3)|p |»l 



(10) 



P± = 



aa 



^(3 



1 - (awpy 

the £ (2) = 0;when |p | > 1 and a > 1 



± 



a 2 -l 



1 - (a( 2 )/?) 2 ^ 1 - « (2) /3 



1/2 



(11) 



;> P- > |Po| 



(12) 



Note that it is impossible to obtain _E = and p e = simultaneously using only the 
second (electron, n b > 0, or positron, nj, > 0) bunch. 

The second root p + will not be reached in the considered cases. 

The electric field E^ has a maximum at 



Pe 



[1 - (aW[3) 2 ] 1/2 

In the considered cases p™ ax is choosen according the condition 



(13) 



\P-\ > |Po| > |p e 
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The length of the second bunch is choosen in such a way that at the end of it E^ = 
-Eq, and p e = p (see Fig. 6). 

The third bunch then can be similar to the first one in order to dump the electric field 
and plasma electron momenta to their initial zero values, in order to dump wake fields 
completely. For this end the length must be found, integrating the equation 

±[(5{l + p 2 e ) - Pe }dp e 



dz 



(14) 



[2(no - 4 2) )] 1/2 (1 + P 2 e ) 1/2 [a - (1 + p 2 ) 1/2 - aW(3p e }V 2 
by two steps-from p to p_ and from p_ to p , taking sign + when p e increases with the z 

and sign — , when p e decreases with the z. In both cases \p e \ > po ^> 1 and the integration 

simplifies essentially, when (1 + p 2 ) 1 / 2 replaced by |p e |;the result of integration, using ([10 

a is 



4(l + /3)|p | 1/2 



(n 



(i) 



n. 



(2)> 



(no 



n 



(2)^ 



a 



P)/3) 



- 1 



1/2 



(15) 



[2(n -n[ 2) )] 1 /2(i_( a (2) /? ))i/2 

In the case, when the second bunch is consisted from positive particles (positrons) 
the length of the second bunch, obtained by the similar calculations, is 



4(l+/3)|p | 1/2 



[2(n + n[ 2) )] 1 /2(i + a(2)^)i/2 



n 



{1) +n { h 2) 



- 1 



1/2 



(16) 



where 



, -(2) 



In the positron case njf* can be choosen as n 



11, 



(i) 



^> n for the effective inversion 



of the electric field E . Comparison of the eqs. (|15"1) and (|l6f) for 

n ,/3 = 1, 



shows, that in this case 



— (2) 

i.e. the dense n b = Uq positron bunch-invertor would be much shorter than electron 
bunch with < l/2n . The length of the first (and third) electron bunch can be 



d% 



n 



2n 



(i) 



< 1 



;i7) 



obtained similiarly or from the results of work [12| for the case, when 



n 



(i) 



1 < -2-<2y,|po| > 
n 
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and is 

di = /3|po| 1/2 = 4 (18) 

The possible largest plasma electron momenta inside the first bunch in the considered 
case differs from that found in by wake wave breaking limit and can be estimated 
from momentum conservation (see also ||10||), 

niW 1} > I"" Pe{z)n e (z)dz = (19) 



4 

For the estimate in (|i~9l) the expression for n e , obtained from continuity equation 

n f3{l + Pe ) l l 2 

is used, which approximated by linear one in interval from n e = no at p e = up to the 
value n e = ^4 for large p e . 
From (HI) 



An 



(i) 

6 



1 < |po| < —^-T, (20) 
n 



and corresponding electric field from @ and (^) is 

|£„|=2( n « 1, ) 1/2 |p„|" 2 <^7" 2 (21) 







As it was shown in |10| the change of the bunch electron momenta Ap, in the first 
approximation, which can be found from the second order equations, is 

Ap = Pbl = -E{z)t (22) 

where E(z) is the electric field in the zero order approximation inside the rigid bunch 
and r = tu p t < 7 1 / 2 . The limit on acceleration time interval follows from the domain of 



validity of the applied multiple scales approach. 

The length of the plasma column needed for the change of the momenta given by (p2j) 
is, in ordinary units 

/ = ct < = £y/2 (23 ) 

Up ZTT 

The acceleration rate (gradient), in ordinary units is 

G = £ iT < 24) 
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4 POSSIBILITIES OF THE EXPERIMENTAL TESTS 

The formulae obtained in sec. 3 for the parameters of proposed device in ordinary units 
and n = n are from fl20|): 



\po\ < ^jmc, (25) 
n 



from (M 



E < 4 7 1/2 £ , 26 

n 



from (HD, (HI) 



4n (1) 

Ap 6 = eE t < 4-^mc7 = — b —p b ; (27) 
n n 

^ < 4^; t <.;V' 2 
\ n 



C = cA Pb - a p - . j f n t I mcr^g (28) 

el et \ rig J e \ J 



The accelerating gradient from (|24[) , (|27D is 

G- % (C) 7 1/2 MV7cm, 
A = 2ss = 6.08 .10- 4 f^) V2 cm 

The acceleration length (^2|): 

l = ct < *P 7 V2. 

~ 2vr ' 



The lengths of the first and third electron bunches (|T8|) , (|25|) are: 



x / (i)\ 1/2 
di = 4 < ^7 1/2 ^- (29) 

The length of the second electron bunch (|15|), (|25| ) is 

^< ^7 1/2 (^) (30) 

( 2 ) ^ ^ (1) 
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The length of the second positron bunch (|T6|), ( p5|) is: 

d p 2 < ^ 7 1/2 (31) 

-(2) (1) ^ 

According to fl26|) it is possible to obtain essential acceleration fields when the previ- 
ously accelerated dense bunches of negative charged particles (electrons) are used. The 
lengths of such a bunches, according to (|29|) , imcreases as ~ 7 1 / 2 , so it seems that they 
must be the bunches from high current accelerators with the energies up to tens of Mev. 
Even at these energies maximal acceleration gradient is high enough, it can exceed the 
acceleration gradient of the ordinary linacs by a several orders of magnitude. 

Consider the experimental possibilities of proof-of-principle experiment on existing 
accelerators or test facilities. It seems that Argonne Wakefield Accelerators (AWA), 
which accelerates very dense bunches, is suitable for above mentioned proposes. 

AWA [|17| is 20MeV electron accelerator which accelerate the bunches of 20p sec full 
width (the length d\ exp = 0.6cm), diameter about 0.1cm [IS|, the total charge of the 
bunch lOOnC, i.e. n^ 1 = 2 ■ 10 14 cm~ 3 . If n is choosen as n = hn^ = 4 ■ 10 13 cm~ 3 , then 
the optimal length of the first bunch from ( |29"D is d\ = 2.4cm, i.e. 4 times larger than the 
experimentally existing one d\ exp = 0.6cm. 

The role of the invertor can serve the plasma column itself. Witness bunch, which 
has lOps duration, less than InC total charge, and energy 4Mev can be used as accel- 
erated bunch, being placed somewhere between first and third bunches. Third bunch 
(damper) has similiar parameters as the first one and must be placed at distance di = 
8 ( J di p& 43cm from the end of the first bunch. Experimentally it is possible 
to put it on the distance d 2 exp = 46cm, because the frequency of the RF-cavities is 
v = 1.3GHz (d2exp = nc/ v = n ■ 23cm, n = 1, 2, 3, • • •). If the length of the first bunch 
could be made 4 times large than existing one and total charge of the bunch will be also 4 
times larger-400nC, then the change of the energy of the witness bunch will be, according 



to 

n (1) 

A£ w < 4-^g 1 = 400MeV > S w = 4MeV. 
n 
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For the existing experimental parameters of AWA:g?i = 0, 5cm, d 2 = 23cm, total 
charge of bunch 70nC, = 1.4 • 10 14 cm results of the numerical calculations, along 
the lines presented in Sec. 2, are shown on figures 8a and 8b. The role of bunch-invertor 
plays the plasma column, with density n as a free parameter, which was been chosen 
from the necessity to invert the field, produced by the first bunch, in such a way, that 
the succesive bunch, with the same parameters as the first one, moving at fixed (by 
RF-frequency) distance from the first one, d 2 = 23cm, will be able to dump the field to 
zero. The sought density of plasma n was found equal to n m 1.17 • 10 13 cm~ 3 , \ p = 
0.973cm, n b /n Q = 11.95. Then the maximum value of the accelerated electric field is 
equal to E max w 100Mv/cm, and E < on the rear part of the second (third) bunch. 
So about the half of the electrons of the second (third) bunch, as it can be seen from 
Fig. 8a, can be accelerated during the time t ~ u;"^ 1 / 2 up to energy Sj 70MeV. 




M.1714X10V, d,=d s =0.5cm, d2=23cm, f=40, l,=0.973cm. Fig 8b phase-plane portrait lor proposed experiment on AWA. 



Proposed experiment can serve as a prof-of-principle one;plasma remains practically un- 
perturbed after the passage of the first and the second bunches, so it is possible to inject 
the next couple of the bunches and also obtain acceleration, again leaving the plasma 
behind the bunches in practically unperturbed state and so on. As it was above men- 
tioned in order to obtain the optimal experimental conditions, it is necessary to enlarge 
the micropulse duration, leaving its charge density the same. This will rise the maximum 
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obtainable electric field and energy gain. Of course, the obtained numbers for E max and 
£ are at some extent tentative;it is neccessary to remember that they are obtained in the 
model of the bunch with the infinite transverse dimensions and uniform charge distribu- 
tion. More detailed calculations must be performed in the course of preparation of the 
experiment. 

However, taking into account the results of the numerical calculation for laser driven 
wake fields in generalized vortex-free plasma |1| , it is possible to adopt the statement of 
authors |Tj|, that "even when driver's width approaches c/u p the 1-D (nonlinear theory) 
predictions are still a good guide for determining the accelerating field strength" . In the 
considered case of AWA the diameter of the bunch is D = 0.1cm < — = tt = 0.15cm, 
so for the existing parameters of AWA predictions of 1-D nonlinear theory can be valid. 

For the improved parameters of AWA, as it was shown, n = 4 • 10 13 cm -3 and D = 
0.1cm > — = ^ = 0.08cm and the situation is more favourable. 

For Accelerator Test Facility (ATF) at BNL the situation in the considered respect is 
worse, due to the very narrow bunches r = 0.03cm and smaller bunch charge Q = 2nC 
p0|| . At any case, the considered problem needs the development of 2-D nonlinear theory 
of wake field generation by the electron (positron) bunches with the finite transverse 
dimension. 

Another possibility to check the predicted effect is based on the use of induction linacs. 
Consider the example of the accelerator with the parameters, which are within the range 
of the possibilities of the existing or planned induction linacs (see e.g. [pl|, [I6[): Energy 
£ = 16MeV 

Pulse current I = 5kA 
Pulse duration r = 2.5nsec 
Diameter of the bunch D = 2cm 

Then the charge density of the first and third bunches will be = 1.6 ■ 10 11 cm~ 3 , the 
plasma density can be chosen as n = |r4 an d then, according to (29, d\ ~ 75cm, the 
distance between second and third bunches, according to ([3(]), is d>2 ~ 13.4m (the plasma 
column plays the role of invertor). The maximum energy gain of the witness bunch of 
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much smaller current (e.g. 0.1 I and r < 2.5nsec), placed in the proper position between 
the first and third bunches, according to fl27D , is 



AS max = 4 ■ 5 ■ 16 = 320MeV. 

The condition of the applicability of the 1-D nonlinear thery in this case is practically 
fulfiled: 

D = 2cm < ^ = 2.3cm. 
~ 2vr 

If the positron bunch can be used as a second invertor bunch with the same current 



I as the first one, then the length ofo, according to fl3~l"|) , can be essentially reduced- 
g?2 = 134.1cm, and the part of positrons from this bunch can be accelerated, up to the 
energies ~ 300MeV. 

It is evident, that in all cases considered, the length of the plasma column must be 
2-3 times larger than the total length of the combined bunch, due to transient effects. 

The increase of the lengths of the bunches with the energy may be considered as tol- 
erable one, when the extended astrophysical objects, as e.g. supernovae, are considered. 
Then proposed mechanism can be applied, for example, even for the explanation of the 
origin of the superhigh energy cosmic rays. This problem (see e.g. ||21|| ) is now the most 
important one for the astrophysics of cosmic rays with the energies £ > 10 17 eV, when all 



known mechanisms of cosmic rays acceleration are uneffective [22 . 

In order to illustrate this possibility consider the following numerical example. Sup- 
pose that the flux of electrons with the initial energy £ \ = 5 • 10 16 eV already ob- 
tained by some effective and commonly adopted mechanism. The flux with the den- 
sity rij = 1.5 • 10 14 cm~ 3 passes through the electron-proton plasma with the density 



uq = 3 • 10 11 cm -3 . Then it followed by positron (|2|) and electron (|3]) fluxes with the same 
densities n& = = Then 
be accelerated up to final energy 



densities n& = = . Then some part of positron and third-electron bunches can 



£f ax j | ^- | s, :m 2 V\ 



n 
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i.e. reaches the atmost limit of the observed cosmic rays. The length of the electron fluxes 
is 140km, the length of the positron bunch is 25km. (Let us mention that up to now 
observed events at such energies are very rare). Considered example is only illustrative 
and aplicability of the proposed mechanism to the problem of the origin of the superhigh 
energy cosmic rays needs, of course, more careful consideration. 

In addition, it is necessary to mention, that the expression for momentum change 
( p7| ) Ap = —eEt, obtained for the first approximation of the developed multiple scales 
method, is valid up to fifth approximation: 

Pb 5 6 -1/2 /oo\ 

— = p b = p b0 + ep bl + e p b5 + e p b6 H , e = 7 , (32) 

mc 

as it follows from the results of 0] . Hence it is possible to assume, that eq. (^7j) is valid 



for the time intervals larger than cu p t < e 1 = 7 1 / 2 used in the present work for estimites. 
But in order to proof this assumption it is necessary to develop more complicated version 



of multiple scales method, than used in [10|, where only three scales are introduced: 



z = z — p Q t, 7~i = et, £1 = ez. 

Enlargement of the validity domain on t needs consideration of the contributions from 
the next scales 

r 2 = e 2 t, t 3 = e 3 t, r 4 = e 4 t 
to the momentum development (P^). 
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